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Abstract 

Two-dimensional (2D) topological insulator (TI) have been recognized as a new class of quantum 
state of matter. They are distinguished from normal 2D insulators with their nontrivial band- 
structure topology identified by the number as protected by time-reversal symmetry (TRS). 
2D TIs have intriguing spin-velocity locked conducting edge states and insulating properties in the 
bulk. In the edge states, the electrons with opposite spins propagate in opposite directions and the 
backscattering is fully prohibited when the TRS is conserved. This leads to quantized dissipation¬ 
less “two-lane highway" for charge and spin transportation and promises potential applications. Up 
to now, only very few 2D systems have been discovered to possess this property. The lack of suitable 
material obstructs the further study and application. Here, by using first-principles calculations, 
we propose that the functionalized MXene with oxygen, M 2 CO 2 (M=W, Mo and Cr), are 2D TIs 
with the largest gap of 0.194 eV in W case. They are dynamically stable and natively antioxi¬ 
dant. Most importantly, they are very likely to be easily synthesized by recent developed selective 
chemical etching of transition-metal carbides (MAX phase). This will pave the way to tremendous 
applications of 2D TIs, such as “ideal" conducting wire, multifunctional spintronic device, and the 
realization of topological superconductivity and Majorana modes for quantum computing. 
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I. INTRODUCTION 


The field of topological insulator (TI) started from the theoretical proposal of two- 
dimensional (2D) TI state in graphene.^®^ Though the band gap of graphene is too tiny 
to be observed,® the conceptual achievement in the band topology has opened the door to 
the field of topological quantum states (TQSs).®®The theoretical proposaPand experimen¬ 
tal verificatiorP of 2D TI in quantum well of HgTe/CdTe have boosted the quick rising of 
the field of TIs. The idea of band topology has been extended to 3D system?®!^ and other 
symmetry protected TQSs.^^ Recently, the band topology in metals, including the topolog¬ 
ical Dirac semimetal,II®IZI Weyl semimetaP^® and Node-line semimetal,l23H23 has also been 
intensively studied. Many of the material realization of these TQSs are firstly predicted 
by theoretical calculations and then confirmed by experimental observations.^^® The bulk¬ 
boundary correspondence of the topological matters is well known now and it is one of the 
most unique properties of them. For example, 2D TI is expected to host quantum spin Hall 
effect (QSHE) with ID helical edge states, namely the electrons in such edge states have op¬ 
posite velocities in opposite spin channels. Thus, the backscattering is prohibited as long as 
the perturbation does not break the time-reversal symmetry (TRS). Such helical edge states 
are expected to serve as “two-lane highway" for dissipationless electron transport, which 
promises great potential application in low-power and multi-functional spintronic devices. 
Large band gap 2D TI is also crucial to realize the long-sought-for topological superconduc¬ 
tivity and Majorana modes through proximity effect.*^® In this point of view, 2D TI is more 
preferred than 3D one, where the backscattering in the surface states is not fully prohibited. 

Compared with the number of well characterized 3D TI materials, fewer 2D TIs have 
been experimentally discovered.!^® The quantum wells of HgTe/CdT^ and InAs/GaStl^ 
are among the well-known experimentally confirmed 2D TIs. Both of them require precisely 
controlled MBE growth and operate at ultra-low temperature. These experimental condi¬ 
tions make further studies hard and reduces the possible applications. There have been 
many efforts to find “good" 2D TIs, which are expected to have the following advantages: 
(1) being easy to be prepared; (2) having large bulk band gap to be operated under room 
temperature or higher; (3) being chemically stable upon exposure to air; (4) being composed 
of cheap and nontoxic elements. The theoretical attempts for predicting good TIs can be 
roughly classified into two categories. 1) tuning the strength of spin-orbit coupling (SOC), 
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i.e., the band gap, based on graphene-like honeycomb lattice, snch as the low-bnckled sil- 
icend^, chemically decorated single layer honeycomb lattice of Sn,I^Gd^and Bi or Sb,l^and 
bncked sqnare lattice BiF®. 2) examining new 2D systems, which might be exfoliated from 
the 3D layered strnctnral materials, snch as ZrTes, HfTe^ and Bi 4 Br 4 .^ Transition-metal 
dichalcogenide (TMD) in IT® and sqnare-octagon haeckelit^®® strnctnre also belong to 
the later category. None of the above has been conhrmed by experiments yet, thongh ZrTes, 
HfTes and Bi 4 Br 4 seem to be very promising, since they do exist experimentally and and 
their single layers are TIs withont any additional tnning. 

Regarding oxide materials, non of them is known to be TI in the experiment, thongh 
there are several theoretical proposals available in the literatnres, e. g., 2D TI in single 
layer of iridate Na 2 lr 03 ,® topological Mott insnlatoi® and Weyl semimetal in pychorelcore 
^ 2 ^ 207 ,® axion insulator in spinnel Osmat^^ and 3D TI in perovskite of YBiOg® and 
heavily doped BaBiOa.® It is generally believed that the strong electronegativity of oxygen 
leads to full ionization of cations and results in ionic bonds with large band gap, which 
makes band inversion difficult. However, the noticeable advantages of oxygen compounds, 
i.e., naturally antioxidant and stable upon exposure to air, have stimulated continuous efforts 
in searching for new oxide TIs. 

In this paper, by using first-principles calculations we demonstrate that the functional¬ 
ized MXenepS® with oxygen, M 2 CO 2 with M=W, Mo and Cr, are 2D TIs. Our phonon 
calculations indicate that the crystal structures are dynamically stable. The band inver¬ 
sion, which is crucial to the nontrivial band topology, is found to occur among the bonding 
and anti-bonding states of M d-orbitals. The results are robust against the use of different 
exchange-correlation functional approximations. The bulk band gap of W 2 CO 2 is as large 
as 0.194 eV within generalized gradient approximation (GGA) and is enhanced to 0.472 eV 
within hybrid functional (HSEOO).®^ Its Z 2 invariant is 1 and has conducting helical edge 
states. Recently, the 2D material MXene have been successfully obtained by the selective 
chemical etching of MAX phases — M„+iAX„ (’^=1, 2, 3 , ...), where M, A and X are a 
transition metal, an element of group 12-14, G or N, respectively.® The bare surfaces of 
MXene sheet are chemically active and are usually terminated by some atoms or chemical 
groups depending on the synthesis process, which are usually fluorine (F), oxygen (O), or hy¬ 
droxyl (OH).®® Therefore, we believe that our proposed M 2 GO 2 can be probably realized 
experimentally in the future and thus, will advance the application of TI greatly. 
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TABLE I: The total energies (in eV per nnit cell) for the six possible confignrations of W 2 CO 2 , 
M 02 CO 2 and Cr2C02- For each confignration, the strnctnre is fnlly relaxed. In each nnit cell, two 
oxygen atoms are reqnired for fnll snrface satnrations. The symbols T, A, and B indicate three 
different absorption sites of oxygen atoms as depicted in Fig. All calcnlations here are performed 
withont spin polarization and within GGA. 


sites of Oxygen 

W 2 GO 2 

M 02 GO 2 

Gr2G02 

TT 

3.162 

2.779 

1.790 

AA 

1.189 

1.112 

0.618 

TB 

1.828 

1.782 

1.370 

TA 

2.024 

1.733 

1.034 

BB 

0.0 

0.0 

0.0 

BA 

0.878 

0.698 

0.329 


II. COMPUTATIONAL DETAILS 

First-principles calculations were carried out by using the Vienna ab initio simulation 
package (VASP)^^*^. Exchange-correlation potential was treated within the generalized 
gradient approximation (GGA) of Perdew-Burke-Ernzerhof type.^ SOG was taken into ac¬ 
count self-consistently. The cut-off energy for plane wave expansion is 500 eV and the k-point 
sampling grid in the self-consistent process was 12x12x1. The crystal structures have been 
fully relaxed until the residual forces on each atom becoms less than 0.001 eV/A. A vac¬ 
uum of 20 A between layers was considered in order to minimize the interactions between 
the layer with its periodic images. PHONOPY has been employed to calculate the phonon 
dispersiorP^. Gonsidering the possible underestimation of band gap within GGA, non-local 
Heyd-Scuseria-Ernzerhof (HSE06) hybrid functionaP^^ is further supplemented to check the 
band topology. To explore the edge states, we apply the Green’s function methocP based 
on the tight-binding model with the maximally localized Wannier functions (MLWE ^^^F^ lgf 
d orbitals of M and p orbitals of G and O as basis set. MLWF are generated by using the 
software package OpenMX! ^^ * ^^! 
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FIG. 1: (Color online) (a) Top view and (b) side view of optimized crystal structure of W 2 CO 2 
and its 2D Brillouin zone. There are three possible sites for oxygen decoration on either side of the 
surface, namely B (on top of C site), T (on top of W in top surface), A (on top of W in bottom 
surface), (c) The triangular prism formed by O and C ions surrounds W. (d) The phonon spectrum 
for optimized W 2 CO 2 . 

III. RESULTS AND DISCUSSION 

The crystal structure of oxygen functionalized MXene is shown in Fig. 1. Bare MXene 
M 2 C is basically a three-layer structure with trigonal lattice. The C atoms form a layer, 
which is sandwiched between two M layers. The in-plane sites of M atoms are (1/3, 2/3) and 
(2/3, 1/3), respectively, of the trigonal lattice of C atoms. The bare surfaces of MXene sheet 
are basically terminated by M atoms and they are chemically reactive. Usually, the surfaces 
are terminated by F, O or OH depending on the synthesis process. This brings a chance to 
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tune the electronic properties of MXene by appropriate surface functionalization.*^^*^ It has 
previously been shown that the MXenes with the full surface functionalizations are ther¬ 
modynamically more favorable than the partial functionalizations,^ where the full surface 
functionalization requires two chemical groups per cell. As shown in Fig. 1, the oxygen 
atom might occupy three different sites on the surface, namely A, B and T on each surface. 
Therefore, there are totally six combinations for decoration of two surfaces as listed in Ta¬ 
ble. Hj For each case, the crystal structure is fully relaxed and the results of total energy 
calculations are summarized in Table. It is observed that the MXeen with BB-type oxygen 
functionalization obtain the lowest energy. In this conhguration the M atom sits inside of 
a trigonal prism formed by the surrounding C and O, which is very similar to Mo atom in 
IH structure of M 0 S 2 . The phonon spectrum of the energetically stable crystal structure is 
calculated and shown in Fig. [^d). Obviously, there is no imaginary frequency, which means 
such structures are also dynamically stable. 

The electronic band structure of W 2 CO 2 is shown in Fig. It is found that there is 
a degenerate band touch point on the Fermi level when SOC is not considered. These 
degenerate states are mostly composed of and dxy orbitals of W atoms as shown 

by the fat-band analysis. If SOC is further considered, it becomes an insulator. Since it 
has inversion symmetry, the parity conhguration of occupied bands at four time-reversal 
invariant momenta in the 2D BZ can be easily obtained and the Z 2 invariant is found to be 
1. This indicates that W 2 CO 2 is a 2D topological insulator with an indirect band gap as 
large as 0.194 eV. The topologically protected conducting edge states have Dirac cone like 
dispersion and connect the bulk valence and conduction bands. 

The total and projected partial density of states, as well as the fat-band plot, clearly 
show that around Fermi level the W d orbitals are dominant around the Fermi level. The 
p orbitals of C and O are mainly located -2 eV below the Fermi level. The hve d orbitals 
of W atom are within the triangular prism crystal held of symmetry. The d^z and dyz 
orbitals are double degenerate and are higher in energy due to the strong hybridization with 
the p orbitals of C and O. The d^-^-y^ and dxy orbitals are also degenerate and they are 
around the Fermi level. The d^'^. orbital has quite weak hybridization with ligand elements 
since it points to the center of the ligand triangle. It is also around the Fermi level and 
slightly lower than dx'^-y^+dxy Therefore, in order to uncover the low-energy physics of 
W 2 CO 2 , we just need to take into account three three d orbitals, dx 2 -y 2 , dxy and ^^ 2 , of the 
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FIG. 2: (Color online) Band structure for W 2 CO 2 calculated (a) without and (b) with spin-orbit 
coupling (SOC). The fat-band are scaled with the projected weight of different atomic orbitals 
within the eigenstates as shown in (a). The comparison with the bands from hybrid functional 
(HSE06) calculation including SOC is shown in (b). Total and projected partial density of states 
are shown in (c). The edge states along lattice constant a is shown in (d). 

W atoms. It is noticeable that there are two W atoms in one unit cell. As shown in Fig. 
the above selected d atomic orbitals form bonding and anti-bonding states. The bonding 
and anti-bonding states of orbitals are even {Aig) and odd (^ 2 ^), respectively. Similarly, 
those from dx 2 _y 2 +dxy are even {Eg) and odd {Eu) but having double degeneracy. The band 
inversion happens between the double degenerate Eg and non-degenerate A 2 U states,which 
brings the nontrivial topology of bands. When the banding effect is considered from T to 
M, the double degenerate Eg states are split and there is no other band inversion happens. 
Further, the SOC introduces the spin degree of freedom and it opens a gap around Fermi 
level as shown in Fig. Due to the heavy W element, the band gap is found to be as large 











































FIG. 3: (Color online) The band inversion mechanism in W2C02. (a) W 5d orbitals are split (b) 
under crystal field (CF) with double degenerate d^ 2 _y 2 +dxy and single degenerate d ^2 orbitals 
around Fermi level E^?. (c) The dimerization of two W atoms in each unit cell leads to bonding and 
anti-bonding states and band inversion between states with different parities (labeled by + and -). 
(d) The band dispersion along F-M does not induce other band inversion, (e) Including SOC opens 
band gap and introduces spin degeneracy in each band, with dashed and solid line representing 
opposite spin channels. 

as 0.194 eV. Considering the possible underestimation of the band gaps within GGA, the 
hybrid functional (HSE06) calculation is used to check whether the above band inversion 
around F is robust. It is found that the band inversion is kept and the band gap is enhanced 
to as large as 0.472 eV. Such large band gap 2D TI has an advantage in observing quantum 
spin Hall effect at room-temperature or higher, which is appropriate for device applications. 

For M 02 GO 2 and Gr2G02, both of them have the similar crystal structure as W 2 GO 2 and 
are also dynamically stable as seen from the phonon spectra shown in Fig. Since Mo and 
Gr have weaker SOG than W, they are expected to have narrower band gap than W 2 GO 2 . In 
fact, both of them show band structure of semimetals by having compensated electron and 
hole Fermi pockets. However, at each k point, one can still hnd well dehned band gap with a 
curved Fermi level. The assumed occupied bands, as denoted with solid lines in Fig. have 
the same Z 2 number as W 2 GO 2 . Therefore, both of them share the same band topology, as 
well as the underlying physics, with W 2 GO 2 . 

Furthermore, we have investigated the correlation effect in d electrons of transition-metal 
M. For quite delocalized 5d electrons in W and Ad in Mo, detailed GGA+G [U is the 
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MojCOj 


CrjCOj 



FIG. 4: (Color online) The phonon spectrum (upper panel) and band structure with SOC (lower 
panel) for optimized stable M 02 CO 2 (left panel) and Cr2C02 (right panel) The assumed occupied 
(unoccupied) bands are denoted with solid (dashed) lines. 

parameter for onsite Coulomb interaction) calculations {U <4.0 eV) show that the correlation 
effect is negligible and the ground state of W 2 CO 2 and M 2 CO 2 are always nonmagnetic. The 
results discussed above are robust. However, for 3d transition-metal Cr, simple GGA+f/ 
{U >2.0 eV) calculation for BB conhguration gives out solution with magnetic properties. 
The above non-spin polarization calculations for Cr2C02 are helpful to understanding the 
band gap dependence on the atomic number in that column, but not realistic for Cr 2 C 02 
itself. However, this gives us the spin degree of freedom in material design based on MXene, 
which is very crucial to seeking Chern insulators hosting quantum anomalous Hall effect 
The detailed studies on the dependence of U value and magnetic conhgurations are left for 
future work. 
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IV. CONCLUSIONS 


Based on the first-principles calculations, we have predicted that a family of oxygen 
functionalized MXene M 2 CO 2 (M=W, Mo and Cr) are 2D TIs. The representative W 2 CO 2 
has robust band inversion, nontrivial Z 2 invariant and large band gap of 0.194 (0.472) eV 
within GGA (HSE06). It might satisfy all the four criteria of a "good" 2D TI. (1) easier 
production process by selective chemical etching method; (2) hosting QSHE at ambient 
condition; (3) high stability and antioxidant upon exposure to air; and (4) low production 
cost and consisting of environmental friendly elements. 

Inspired by these Endings, one can naturally think about obtaining more 2D TI candidates 
in other functionalized MXenes. The large number of MAX (more than brings many 

possibilities of MXen^SQl ^nd huge space in finding topologically nontrivial materials. The 
possible changes would include replacing with F“ or varying transition- 

metal M, replacing G with N or tuning the number of layers n in MXene M„_|_iG„, etc. 
The material design or property tailoring with a single or any combination of the above 
changes will lead to more and better 2D TIs or Ghern Insulators. 
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